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ABSTRACT 

In this paper we discuss the recently discovered radio transient in the nuclear region 
of M82. It has been suggested that this source is an X-ray binary, which, given the 
radio flux density, would require an X-ray luminosity, Lx '^ 6 x 10^^ erg s~^ if it were a 
stellar mass black hole that followed established empirical relations for X-ray binaries. 
The source is not detected in the analysis of the X-ray archival data. Using a 99 % 
confidence level upper limit we find that Lx ^ 1-8 x lO'^^ergs"^ and 1.5 x lO'^^ergs"^, 
using powerlaw and disk blackbody models respectively. The source is thus unlikely 
to be a traditional microquasar, but could be a system similar to SS433, a Galactic 
microquasar with a high ratio of radio to X-ray luminosity. 
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1 INTRODUCTION 

M 82 is prolific star-for ming galaxy at a distance of 3.6 Mpc 
l|Freednian et al.l 1 1994 ). Its proximity therefore makes it 
an ideal place to study star formation. To this end, M 82 
has been well monitored by vario us observatories at sev- 
eral wavelengths for d ecades (e. g. 'TcIcsco & Har pej 198(3: 
Kronberg et al] Il98ll: iBland fc T uUv 1988; Muxl ow et al 



1994 



Strickland et all 1 19971 : IWeifi et al.ll2001l : IFenech et al 



20081 ). In particular, radio observations have revealed ap- 



proxima tely 60 compact radio sources in the central region 
of M82 iMcDonald et al.|[20o3 ). A quarter of these sources 
are of unknown origin. In addition to compact sources, fre- 
quent radio observations have als o found transient sources of 
an undetermined na ture in M 82 (jKronberg &: Sramekll 19851 : 
IMuxIow et al.l 1 19941) and monitored the evolution of radio 
supernovae (e.g. iMuxlow et all 1 1994 iBeswick et all l2006l : 
iFenech et aLlbOOsf T 

The next generation of radio telescopes should reveal 
many more such transient sources. These observatories will 
have the ability to detect a 0.1 mjy source within minutes. 
We will then more easily be able to observe transient pop- 
ulations not only in star-forming galaxies like M 82 where 
these systems are abundant, but also in galaxies with low 
star-formation rates like ellipticals. By attempting to deter- 
mine the nature of the unknown transients such as those in 
M 82, we will make the task of studying the large number 
of transients expected to be observed in the future much 
easier. 

In this paper we investigate further the most recently 
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discovered radio transient of unknown origin in M82 and 
show that this source is unlikely to be a normal microquasar. 
We discuss the possibility that the M 82 source could be an 
extragalactic analogue of S S 433, a very unusual Galactic 
microquasar (|Margonl 1 19841 . and references therein) . Both 
sources have low X-ray and high radio luminosities. In ad- 
dition, SS433 possesses precessing, helical jets. SS433 also 
has highly blu e- and redshifted op tical emission lines asso- 
ciated with it (JMargon et al.lll979l ). If the M82 transient is 
an extragalactic SS433 like microquasar, it will be only the 
second such sourc e, after the microq uasar S 26 in NGC 7793, 
to be discovered (|Pakull et al.ll201ol ). 



2 PROPERTIES OF THE SOURCE 

The radio source was first discovered with the Multi-Element 
Radio-Linked Interferometer Netwo rk (MERLIN) early in 
May 2009 by iMuxlow etlH (|2009l ) as part of a continu- 
ous monitoring o f the r ecent supernova SN 2008iz in M82. 
iBrunthaler et al.l l| 20091 ) later also reported a detection of 
the source on 2009 April 30 with the Very Long Baseline 
Array (VLB A). 

For the last 15 to 20 years M82 has been observed at 
radio frequencies at intervals of six months to one year. In 
that time, only two other transient so urces of unknown na- 
ture have been detected i n the galaxy IjKronberg fc Sramekl 
1 19851 : IMuxIow et al.lll994l). This most rece nt transient source 
was not detected bv lMuxlow et al.l (|2009|) in the period 2009 
April 24th t o 27th to a 3 g upper limit of < 0.2 mJy/ beam 
at 4.9 GHz. IBrunthaler et all (|2009i ) found that the source 
was not detected in their observations taken on 2009 April 
8 (3cr limit of 0.6 mJy at 22 GHz) and 2009 April 27 (3cr 
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limit of 0.9 mjy at 4 3 GHz and 0.7 mjy at 22 GHz). More- 
over, IMuxIow et al.l l|2010r ) report that the source has not 
been detected in other wavebands previously or con currently 
with their radio detections, including infrared ( Fraser et al.l 
l2009f ). optica l (Matilla, private com munication) and X-ray 
observations iJKong fc Chianell2009l 'l. 

The source is situated in the nuclear region of 
the galaxy, at a position of RA 09''55™52^5083, Dec. 
69°40'45".410 (J2000) with 5 mas error in each coordinate. 
The VLB A observations show flux densities of 1.1 mJy and 
0.5 mJy at 1.6 GHz and 4.8 GHz respectively. The observa- 
tions from the first week of 2009 May revealed a peak radio 
flux of between 0.6 and 0.7 mJy at 4.9 GHz. Follow up obser- 
vations showed tha t the flux remained r oughly constant over 
the next 150 days IjMuxlow et al.ll2010l ). The radio spectral 
index, a, also remained unchanged at —0.7 during this time, 
where F,y oc u°' . MERLIN data taken in 2010 January (over 
240 days after the initial detection) show that the source is 
roughly 15 mas in size. 



3 X-RAY OBSERVATIONS 

The nucleus of M82 was observed with Chandra ACIS-S 
on 2009 April 17 (ObsID 10025) and 2009 April 29 (ObsID 
10026). The data files were downloaded from the Chandra 
on-line archive. The data were checked for background flares 
and spectra were produced using the CIAO script psextract. 
We used a circular region of radius 2" centered on the radio 
source to extract the source and background spectra from 
ObsID 10026 and 10025 respectively. 

The spectra were analysed using the Interactive Spec- 
tral Interpretation System (ISIS), version 1.6.1. As ObsID 
10025 was taken well before the source was flrst seen in ra- 
dio, we used this observation as the background spectrum 
for our analysis. Given the initial, sharp radio variability, we 
expect that the variable component of the X-ray flux should 
be comparable to the X-ray flux itself. We were unable to 
use the standard detection experiment approach because of 
the strong diffuse X-ray emission in the centre of M 82. The 
data were grouped into bins of 50 counts to increase signal- 
to-noise; binning to the standard signal-to-noise ratio of 5 
produced too few bins for fitting. Only data in the range 
0.3— lOkeV were included in the spectral fitting. Due to the 
low number of counts only two spectral models were used, 
namely disk blackb ody (diskBB) and power law (PL) (see 
iMitsuda et"aLlll989l ). We fixed the PL index, F = 2 and the 
diskBB temperature, Ti„ — 1 keV. The Gehrels statistic was 
used to calculate error s as it gives a more reliable fit to data 
with low count rates l|Gehrelall986l ). Both the diskBB and 
PL models give similar results. 

The spectral fitting of the data yielded no clear 
detection of significant X-ray flux. In fact, the back- 
ground spectrum had a higher number of counts than 
the source spectrum. We calculated the upper limits cor- 
responding to a 99 % confldence level for the X-ray lu- 
minosity and found that the unabsorbed luminosities are 
1.8x10^^ ergs"^ and l.SxlO^'' ergs"^ for PL and diskBB 
respectivel y, using Galactic column density, uh = 4 x 
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Figure 1. The PL model fit (99 % confidence level upper limit) 
of X-ray (0.3-10 keV) data for the M82 source. TheAx^ contri- 
butions from each bin are shown in the bottom panel, with the 
sign indicating when the data are above or below the model. 



Fig.m The data are in the range 0.3— lOkeV. The bottom 
panel shows the Ax^ value from each bin. 

The source is situated at a local minimum in fore- 
ground column density with a value of 3. 5x10^^ cm~^ esti- 
mated from CO maps (jWalter et al.ll2002l ). Using this mea- 
sured column density, the upper limits corresponding to a 
99 % confidence level for Lx in the 0.3 — 10 keV range is 
4.9x10^^^ erg s"^ (PL) and 3.1x10^^^ ergs"^ (diskBB). The 
X-ray non-detection is therefore unlikely to be caused by 
absorption in M 82. 



PROPERTIES OF GALACTIC BLACK HOLE 
BINARIES 



-.20 , 



10^" cm"^ JDickev fc Lockman|[l990h . The spectral fit for 
the 99 % confidence level limit of the PL model is shown in 



IMuxIow et al.l (|20ld ) discuss various possibilities regarding 
the nature of the source, including that the source is an ex- 
tragalactic microquasar. Other ideas put forward are that 
the source might be an AGN in the centre of M 82 or an un- 
usual radio supernova. However, they also note that both 
these explanations are unlikely. If the source is an AGN 
it would require the supermassive black hole to be signif- 
icantly displaced from the dynamical centre of the galaxy. 
The source is also too faint to be a Type I supernova and 
the unchanging spectral index argues against young radio 
supernovae in general. 

In this paper we focus on the microquasar scenario. We 
use the properties of black hole binaries (BHBs) to predict 
various parameters, such as X-ray luminosity, for our source. 

An unambiguous correlation has been found between 
X-ray luminosity, Lx, and radio luminosity, Lradio, for stel- 
lar mass black holes in the hard X-ray state (always the 
case below 1% Eddington), of the form Lradio oc Lx where 
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0.6 s: b < 0.7 (|Corbel et al.ll2003l : iGallo et al.ll2003l . l200d ). 
At higher Eddington ratios, X-ray binaries can enter phases 
of radio flaring associ ated with state transitions. As dis- 
cussed in lFender et al.l ((2004) (see for example their Fig. 2), 
in such states the radio/X-ray correlation is still broadly 
consistent with the data, albeit with a larger scatter (of 
around 1 dex in Lradio, corresponding to a scatter in Lx 
of -^ 1.4 dex). This radio/X-ray correlation has been shown 
to extend, with the addition of a mass term, M, to active 
galactic nuclei (AGN) , in a tight corre lation for AGN at 
low Eddington ratios (|Falcke et al.ll2004h . and with a simi- 
lar correlation but broader scatter when including a sample 
of AGN with no restriction on Eddingto n ratio or state, and 
which thus must include flaring sources (JMerloni et al.ll2003l . 
hereafter MHdMOS). The best flt to the MHdMOS 'funda- 
mental plane' is of the form Lradio oc L'^^M^'^, and has a 
scatter of about 2 dex in Lradio- The current state of knowl- 
edge is, therefore, that for 'normal' X-ray binaries and AGN 
there is a common radio/X-ray correlation which is tighter 
in hard states but which follows a similar, but less precise, 
correlation even in flaring states. There are of course addi- 
tional caveats if the radio observations are very sparse (e.g. 
whether or not the flare peak was caught), but in the case 
of the M82 transient the sampling is good. 

Assuming that the M82 source is a normal BH X-ray 
binary transient following the Corbel/Gallo relations, the 
measured radio luminosity implies Lx ~ 6 x 10*^ erg s~^, 
far in excess of our established upper limit (see Fig. [3]). 
Using the MHdMOS fundamental plane and assuming the 
M82 source is Eddington-limited, then Lx = ^Edd = 
1.3 X 10^* (M/M©) erg s"^ implies M = 1.9 x 10^ Mq and 
Lx = 2.5 X 10*^ergs~^. This inferred luminosity suggests 
that the source should be very bright in X-rays and that it 
could be an intermediate mass BH system. For the Edding- 
ton limited case Lx oc M and thus Ljadio ex Lx^, leading 
to a scatter to 1.4 dex in Lx, the same as that estimated 
for the radio/X-ray correlation. This scatter introduces an 
uncertainty of a factor of about 30 into the implied Lx. If 
Lx is reduced by a factor of 30, the resulting luminosity is 
still several orders of magnitude higher than the upper limit 
estimated from the spectral fits. 

It should be noted that Lradio/Lx could sometimes be 
higher for HMXBs than for LMXBs. HMXBs have dense 
stellar winds that might interact with the jets from the BH 
(see e.g. lRomero c t al. 2003). Any jet-wind interaction could 
give rise to synchrotron radiation thereby increasing the ra- 
dio emission of the system. The mass transfer rate, M, of 
the binary system also increases with the mass of the donor 
star (see §5). Additionally, HMXBs in Roche lobe overflow 
may have donor stars sufficiently massive to allow unstable 
mass transfer. If this happens, the accretion rates can be- 
come highly super-Eddington. When M is strongly super- 
Eddington, the accretion disk is expected to become geomet- 
rically thick a nd the X-rays from th e BH are scattered and 
absorbed (e.g. lBegelman et al.ll2006l V Both the low radiative 
efficiency in the X-rays and the high radiative efficiency of 
the jets in radio then serve to increase Lradio/Lx. 

Now we compare the M 82 source to other non-standard 
BHBs to try and establish whether or not our source could 
be such a system. Cyg X-3 is the brightest quasi-persistent 
Galactic BHB at radio wavelengths with a peak fiux den- 
sity of approximately 20 Jy at 2.25 GHz. When scaled to 
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Figure 2. Top panel: A 150 day light cu rve of Cyg X-3 taken 
with the GBI at 2.25 GHz dWaltman et aft 19951. The horizontal 
line indicates the flux density of our source scaled to the distance 
of Cyg X-3. Bottom panel: A light curv e of SS 433 taken with the 
GBI at 2.7 GHz fsee lFiedler et al.|[T987l . for data up to December 
1985 and procedure for subsequent data downloaded from the 
GBI website). 



the distance of Cyg X-3 (rf^Qkpc. IPredehl et ai]l2000l ) the 
peak flux density of our source is 110 Jy, nearly ten times 
that of the brightest Cyg X-3 flare. Thus this source would 
be exceptionally radio bright if it is a microquasar. More- 
over, the spectral index of Cyg X-3 varies greatly during 
flaring episodes, with —0.4 < a < 1.8 (jWaltman et al.l 
1 1995). It is then unlikely that Cyg X-3 and the source in 
M 82 are in the same accretion state. Figure[2] (top panel) 
shows a radio (2.25 GHz) light curve of Cyg X-3 for a pe- 
riod of 150 days taken with the Green Bank Interferometer 
(GBI) (jWaltman et al.lll995l ). The horizontal line indicates 
the peak flux density of our source scaled to the distance of 
Cyg X-3. 

Most other Galactic microquasars have radio flar- 
ing events that show a rapid increase in flux density 
over periods of hours to days; the flares then decay 
back to quiescent levels within a few days to weeks 
(see e.g. iH annikai nen et al.l l20(3ll: iBrocksopp et al.l I2OO2I : 
iFender et al.ll2002l : IStevens et al.ll2003l ). The flux density of 
the M 82 source increased by at least a factor of five over a 
period of a few days. Unlike most other microquasars, how- 
ever, the source then remained at approximately constant 
flux for several mon ths, with flux dens ity va riations of less 
than a factor of two lMuxlow et al.l (see|2010|. Fig. 2). 

SS 433 is an Galactic exotic microquasar that in some 
respe cts exhibits ra dio emission unlike that of other BHBs 
(e.g. iFabrikal l2004l l. A light curve taken with the GBI 
at 2.7 GHz for this system is shown in Fig. [2] (bottom 
panel). The d ata u p to 1985 December were taken from 
iFiedler et al.l l| 19871 ). Subsequent data were downloaded 
from the GB I web site q the procedure is the same as in 
iFiedler et all (|l987l 'l. The fight curve of SS433 is similar to 
that of our source as its flux density stays roughly constant 

^ http://www.gb.nrao.edu/fgdocs/gbi/gbint.html 
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Figure 3. Plot showing radio/X-ray correlation for hard state 
BHBs . The squares indicate the BHB data from ICalvelo et alj 
(|2010|); the X-ray data are in the 1— lOkeV range and the radio 
data are calculated at 5 GHz. The triangle shows a representative 
SS 433 data point (2— lOke V and 4.9 GHz for X-ray and radio data 
respectively; iMargorJI 19841 . and references therein) . The asterisks 
indicate Cyg X-3 data (2—11 keV an d 4.9 GHz for X-ray and radio 
data respectively: iGallo et al.ll2003l) . The filled and open circles 
are the actual and estimated (from the radio/X-ray correlation) 
data for the M82 source respectively, with Lx in the 1 — lOkeV 
range and Lradio calculated at 4.9 GHz. The errors have been 
omitted from the plot for clarity; see MHdM03 for error estimates. 



over a period much greater than the decay time of most mi- 
croqu asar fiares. In addit ion, the spectral index of SS 433 is 
—0.5 (IDubner et al.lll998r . similar to that of the M 82 source. 
S 26 is alre ady thought to be an extragalactic ana- 
logue of SS433 (jSoria et al.ll20ld '). In contrast to the point- 
like M82 transient, S 26 has extended structures of radio 
lobes and X-ray and radio hotspots nearly 300 pc apart en- 
veloped in a cocoon of gas that has been inflated by the 
jets. The radio spectral index varies across the system. In 
the lobes —0.7 ^ a ^ —0.6, it flattens across the cocoon 
with —0.4 ;$ a ^ and is inverted near the base of the 
jets where < a < 0.4. The core of S 26 has an X-ray lu- 
minosity of ~ 7 X 10^^ ergs"^ (0.3 - lOkeV), but no radio 
emission has been detected from the core (3 a upper limit '^ 
3x10'^'' ergs~^). This more complex structure suggests that 
if S 26 and the M 82 source are binary systems with simi- 
lar components (i.e. similar donor star and accretor), then 
S26 has been undergoing mass transfer for a much longer 
duration so that it has had a chance to power strong lobes. 



5 DISCUSSION 

In Fig.|3] we present the observed upper limit on the X- 
ray luminos i ty and the X-ray luminosity predicted from the 
iGallo et al.l (|2003l ) correlation versus the radio luminosity 
on the plot of the observed X-r ay and radio l umino sities 
for BHBs, using the data from ICalvelo et al.l (|2010l ). For 
the purpose of this comparison we recalculate Lx for the 
M82 source to be in the 1 — lOkeV range and Lradio is cal- 
culated from the flux density data at 4.9 GHz. Data for 
SS 433 (lMargonlll984 and references therein) and Cyg X- 
3 (|Gallo et al.l 120031 ) are also included in the plot. The X- 
ray luminosity for the M 82 source is lower than any other 



transient Galactic BHB, yet the radio luminosity is excep- 
tionally high. This, along with behaviour of the Ught curve, 
makes it unlikely that the source is a normal microquasar. 
Like the M 82 source, SS 433 is also relatively faint in X-rays 
compared to radio, with Lx ~ 3 x lO^^ergs"^ (2— lOkeV) 
and Lradio ~ 7 x 10^^ ergs"^ at 4.9 GHz. The M 82 transient 
source could therefore possibly be an extragalactic analogue 
of this system. 

The fact that the M 82 source has a significantly higher 
radio luminosity than SS433 could be explained ,perhaps, 
by a more massive donor star. The mass of t he donor star in 
SS43 3 is thought to be only about 12 M© (JHillwig fc GiesI 
2008). This mass is at the lower limit for systems powere d 
by thermal timescale accretion onto a BH (jKing et al.ll2000l ). 
It is worth bearing in mind that at the distance of M82, 
SS433 and the W50 nebula would be difficult to distinguish 
from a supernova remnant. Additionall y, the flux density o f 
W50 (d = 5.5 kpc) is 70 Jy at 1.4 GHz |Dubner et al.lll998l ). 
At the distance of M82, the flux density would be 160 /xJy 
spread over an angular scale of 7". At these flux densities 
such sources will be impossible to detect as they are well 
below the confusion li mit for M 82 (1.45 mjy with 1" angular 
resolution at 4.9GHz: rKronberg et al.lll98ll) . Thus only the 
brightest and/or most compact SS433-like systems will be 
detectable in galaxies like M82. 

Strong blue- and redshifted optical emission lines like 
those seen in SS 433 could potentially be used a diagnos- 
tic tool to verify the nature of our source. SS433 has an 
R band magnitude of 12 mag; the Ha emission line is 10 
times brighter than the continuum and 5 000 kms^^ wide. 
At the distance of M 82, assuming similar intrinsic luminos- 
ity and extinction, the continuum of our source would not 
be detected with ground based instruments, but the Ha 
line should be. Given that the radio luminosity of the M82 
source is many times greater than that of SS 433, it would 
not be surprising if the M82 source were also more luminous 
optically, increasing the chances of detecting the emission 
lines. However, if the M 82 source suffers from more extinc- 
tion there may also several bright infrared emission lines 
(Paschen a and B rackett 7 to B15) that could be used to 
identify it instead (|Thompson et al.lll979l ). 



6 CONCLUSIONS 

We have shown that the recently discovered transient in 
M82 is relatively faint in X-rays given its radio luminosity. 
The lack of variability in its radio Ught curve also make 
it unlikely that the source is a normal BH transient. The 
source could very well be the extragalactic analogue of 
SS433 or a younger version of S 26. 
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